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ABSTRACT 

We examine the projected correlation of galaxies with mass from small scales (<few hundred kpc) 
where individual dark matter halos dominate, out to 15 Mpc where correlated large-scale structure 
dominates. We investigate these profiles as a function of galaxy luminosity and redshift. Selecting 0.8 
million galaxies in the Deep Lens Survey, we use photometric redshifts and stacked weak gravitational 
lensing shear tomography out to radial scales of 1 degree from the centers of foreground galaxies. We 
detect correlated mass density from multiple halos and large-scale structure at radii larger than the 
virial radius, and find the first observational evidence for growth in the galaxy-mass correlation on 
10 Mpc scales with decreasing redshift and fixed range of luminosity. For a fixed range of redshift, 
we find a scaling of projected halo mass with rest-frame luminosity similar to previous studies at 
lower redshift. We control systematic errors in shape measurement and photometric redshift, enforce 
volume completeness through absolute magnitude cuts, and explore residual sample selection effects 
via simulations. 

Subject headings: cosmology: observations - gravitational lensing - dark matter - large-scale structure 
of Universe, galaxies: evolution 



1. INTRODUCTION 

The presence of dark matter in the Universe is well- 
established and contributes significantly to structures 
ranging from galaxies to superclusters. Less understood 
is the distribution and evolution of dark matter corre- 
lated with galaxies over larger volumes beyond the galaxy 
virial radius. The WMAP cosmic microwave background 
(CMB) result for fi m = 0.27±0.03 (jKomatsu et al.ll2011fl 
is over twice that derived from N-body model fits to lens- 
ing studies of the inner mass profile of lum inous galaxies 
toget her with cosmic luminosity density (jMasaki et al.l 
2012). This is not surprising, since about half the 
dark matter (DM) is expe cted to be in virialized halos 
(|Fukugita fc Peebles 2006). Detailed measurements of 
the mass distribution on large scales around galaxies, or 
"galaxy-mass correlation," as a function of galaxy prop- 
erties can thus be a diagnostic of structure formation 
and evolution. A universal mechani sm for hierarchical 
struct ure formation was developed bv lPress fc Schechterl 
(11974ft. drawing o n bot tom-up s tructu re formation ideas 
of iPeeblesI (fl965h and iSaslawl (fl96l . More generally, 
growth of very large scale mass structures is cosmology 
depend ent, and is one of the p robes of the physics of dark 
energy (|Albrecht et al J 12006ft . 

Weak gravitational lensing (WL) is the only direct 
probe that can both measure the mass profile associ- 
ated with galaxies over a wide range of radii crossing 
the virialization and turnaround scales and does not re- 
quire assumptions about the dynamical state or baryon 
content of the system in question. WL is an inherently 
statistical technique: many galaxies are required. The 
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dimensionless shear signal ranges from less than 0.1% 
(large-scale structure) to 1% (galaxies) to 10% (clus- 
ters of galaxies). The WL signal from foreground DM 
halos lensing more distant galaxies is commonly known 
as "galaxy-galaxy lensing" (GGL) and is measured by 
cross-correlating the positions of foreground DM halos 
(as traced by their resident galaxies) with the lensing- 
induced shear of background galaxies. Most GGL re- 
sults in the literature focus on the inner parts of the DM 
halo mass profile inside the virial radius, which is dom- 
inated by contributions from individual halos. Observa- 
tions of lum inosity scaling of the h alo mass have been 
report e d bv iHoekstra et a~l (12005ft . IMandelbaum et aT 



( 2006) . iKleinheinrich et alT(|2006ft . and Ivan Uitert et al 



(|2011ft . Recent studies have also tackled the prob l em of 
redshift evolution of h alo mass ()Lagattuta et al.l 120101 : 
iLeauthaud et al.l [2012). However, the study of GGL 
outside the virial radius is difficult from both obser- 
vational and theoretical standpoints. In observations, 
residual systematics can masquerade as the le nsing sig- 
nal on large scales (jMandelbaum et"aLll2005aft . On the 
theory side, we need detailed N-body simulations and 
ray-tracing studies or compre hensive statistical methods 
to interpret the signal (e.g. IMandelbaum et al.l l2005bl : 
IHavashi fc Whlt^l200l iMasaki et al.H2012ft . 

With five 4 square degree deep fields imaged in BVRz', 
the Deep Lens Survey offers the possibility of measur- 
ing galaxy-mass correlations tomographically over a wide 
range of projected radii well beyond the halo virial ra- 
dius, and over a useful range of redshift. In this paper, 
we investigate how the average galaxy-mass correlation 
over a wide range of radii varies with luminosity and red- 
shift up to z=0.75. We show the first observational ev- 
idence for growth of projected galaxy-mass correlations 
on 10 Mpc scales. Section 2 briefly describes the formal- 
ism for weak lensing; Section 3 describes the data with 
particular attention to the measurement of both image 
shapes and photometric redshifts; Section 4 describes our 
investigation of photometric redshift error propagation; 
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we present the results in Section 5 and compare them 
with previous observations and simulations; and we dis- 
cuss and summarize our results in Section 6. The Ap- 
pendix contains further details of systematics tests and 
correlation matrices. Throughout, we assume a ACDM 
universe with H — 70 km Mpc s , = 0.27, and 
Q\ = 0.73. Distances are given in units of comoving 
Mpc. 

2. FORMALISM 

2.1. Gravitational Lensing 

The projected mass overdensity of a lens is related to 
the induced shear of source galaxies through 

7T E crit = E(< R) - E(fl) = AE (1) 

where R is the transverse separation between lens and 
source on the sky, and 7t is the tangential shear az- 
imuthally averaged over an annulus with radius R cen- 
tered at the lens. AE is the mean excess projected sur- 
face density, which is defined as the difference between 
the average surface density within radius R, E(< R), 
and the projected sur face density at that radius E(i?) 
(|Miralda-Escuddll991[ ). In this paper, we refer to AE as 
a probe of galaxy-mass correlations, not to be confused 
with the 3D galaxy- mass correlation function £ gm . AE 
and £ gm are related in the following way: 



E(J2)=y PtgmWR? + X 2 )dx 



E = — / E(R')R'dR' 
R Jo 



(2) 



(3) 



p is the mean matter density, and x is the line of sight 
distance. The critical surface density E cr i t is defined as: 



J crit 



AttG D l D ls (1 + z l ) 2 



(4) 



where D^, Dg, and Dls are the angular diameter dis- 
tances to the lens, source, and between lens and source. 
We will notate lens and source redshifts as z^ and zg, 
respectively. The factor of (1 + zl) 2 arises from our use 
of comoving coordinates. This choice is motivated by the 
relatively long baseline in lens redshift and radial sepa- 
ration scales that we investigate in this paper. 

If we correct for the systematic effects introduced in 
the observing process and assume that gravitational lens- 
ing is the only phenomenon systematically distorting the 
shapes of the observed background galaxies around the 
position of foreground galaxies in the weak limit where 
both shear 7 and convergence k = E/E cr it are much 
smaller than 1, 7t is related to the observed ellipticity 
er via: 

e T = 7T + e ^ (5) 

where e^ is the ellipticity component perpendicular to 
the line transversely connecting the lens and source and 
e™* is the intrinsic ellipticity. Note that the definition of 
ellipticity used here is e = (a — b)/(a + b) where a and 
b are the semimajor and semiminor axes, respectively. 
When averaging over an ensemble of randomly oriented 
source galaxies, the (e™ ') term drops out, leaving: 



While Equation [5] is only valid in the limit e™' <C 1 , 
we account for this by normalizing to simulations as de- 
scribed in Section [3~T1 

Our estimator for the differential surface mass den- 
sity for stacked lens galaxies is a lensing signal-to-noise 
variance weighted average summed over lens-so urce pairs 
(jSheldon et al.|[200l iMandelbaum et al.ll2008D : 
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where weights wTi depend on angular diameter distances 
(see Equation |4{, the shape noise ctsn, and the elliptic- 
ity measurement error a e . The shape noise is the rms 
ellipticity per component for the source sample. In prac- 
tice, we measure AE for each lens-source pair and bin the 
values by the projected comoving separation between the 
lens and source. The relationship between the estimated 
differential surface mass density AE and the true quan- 
tity AE is investigated using simulations in Section FTTl 
and for the rest of the paper the tilde is dropped. 

2.2. Halo Mass Profiles 

The mass density profiles of dark matter halos are com- 
monly modeled using Navarro- Frenk- White (NFW) func- 
tions, which were mo tivated by dark matter simulations 
(|Navarro et al.lll996Q . 



where 



Pnfw(^) 



5 = 



Sp c (z) 



(t-)( 1 + 7 1 ) 2 



200c J 



(9) 



3(ln(l 



is a characteristic density, p c (z) is the critical density at 
redshift z, and r s is the mass scale radius. The concen- 
tration c = r2oo/r s , and r2oo is the radius at which the 
mean enclosed mass density of the halo equals 200 x p c . 
M200 is the mass enclosed within r2oo, and we use the 
terms "M200" and "halo mass" (Mhaio) interchangeably, 
and we note the specific cases where M^aio means some- 
thing other than M2oo- 

We obtain the surface mass density profile Enfw(-R) 
by integrating the three-dimensional density over the 
line of sight. A second integration yields the mean sur- 
face mass density inside the projected radius: Enfw(< 

R) . T he expressio ns for these quantities c an b e 

found inlBartelmannl (|1996D : iWright fc Brainerdl (|2000D : 
lYang et all (|2003D . 

The enclosed mass, radial scale, and concentration are 
degenerate. To reduce the degenera cy, we fix the con - 
centration using the relation found in lDuffv et al.1 ()2008f ) 
base d on N-body simulatio ns using the WMAP 5 cosmol- 
ogy (jKomatsu et al.ll2009H . effectively turning the NFW 
model into a single parameter profile. 



c(M 200 ,z) = 5.71( 



Mo 



2 x lO^h^Mp 



{It) = (er) 



(G) 



-0.084 (1 + z) -0.47 

(10) 
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In this work, we do not fit a full halo model, which can 
be used to isolate contributions from individual halos, 
satellite halos, and neighboring halos as is done in, e.g . 
iMandelbaum et all (|2005bt ) and Ivan Uitert et all (j201lj ). 
As it is a useful tool for interpretation purposes, we in- 
vestigate the halo model as a fitting function in future 
work and simply comment on how its use might impact 
our results in this work. 

3. DATA 

The Deep Lens Survey (DLS) consists of ~100 nights 
of BVRz' imaging in five widely separated fields DLS 
F1-F5, each subtend ing 4 square degrees of the sky 
(|Wittman et al.|[2002t ). Each field is composed of nine 
"subhelds," each subtending an area (40' x 40') slightly 
larger than the camera field of view and covered with 
dit hers of ~200". FI and F2 were observed by Mosaic- 
1 (|Muller et all [l99l at the NOAO/KPNO 4-m May- 
all Telescope, and F3-F5 were observed by Mosaic-2 at 
the NOAO/CTIO 4-m Blanco Telescope. The DLS was 
designed to go deep enough at high S/N to yield good 
photometric redshifts and to measure galaxy shapes at 
low surface brightness, while being wide enough to aver- 
age over sample variance and to contain enough galaxies 
to minimize shot noise in precision studies of mass. The 
observing strategy was tuned specifically for lensing stud- 
ies: imaging in the R band only when the seeing was 0.9" 
or better. The effective exposure time in the R band is 
18, 000s, and the effective exposure time in the other three 
bands is 12, 000s. We performed inte rnal photometric cal- 
ibrat i on using the ubercal method (jPadm anabha n et al.1 
120081: IWittman et al.ll20lTI ) and external cali bration us- 
ing lLandoltl (If 9921) standar d stars and SDSS (|York et al.l 
I2000t IWittman et al.|[20lfl ). Details of the photometric 
stacking procedure can be found in the DLS technical 
summary paper (Wittman et al., in prep). We created 
object catalogs using SExtractor in du al-image mode us- 
ing t he R-band image for detections ([Bertin fc Arnoutsl 
1996). All magnitudes in this paper are SExtractor 
MAG_AUTO (henceforth referred to as run) calibrated 
to the Vega system and have been corre cted for galac- 
tic ext inction using the reddening maps of lSchlegel et al.1 
(1998) and for point spread fu nction (PSF) var iation us- 
ing the ColorPro algorithm of iCoe et al.1 (|2006h . To con- 
vert the R band magnitude to the AB system, one must 
add 0.200 The survey is 50% complete for object recov- 
ery to niR = 26, nig y = 25.5, and m z = 24.5 where the 
subscripts B, V, and z refer to the B, V, and z bands, 
respectively. After conservative quality and S /N cuts de- 
scribed below, we utilize a subset of ^800,000 galaxies. 
The analysis processes for shear and photometric red- 
shift measurement are described in Sections 13.11 and 13.21 
below. 

3.1. Shapes 

In order to correct the observed ellipticities to the 
intrinsic ellipticities in Equation [SJ we must account 
for contributi ons from i nstrum ental and atmospheric ef- 
fects. As in lJee efaLl (|2007l ). we use an interpolation 
of principal components characterizing star ellipticities 
and sizes measured on individual exposures to produce 

1 See http://dls.physics.ucdavis.edu/calib/vegaab.html 



a PSF model for the co-add image. To summarize the 
procedure, we select high S/N isolated stars using an 
iterative algorithm based on half-light radius vs magni- 
tude, measure principal components (eigen-PSF) and co- 
efficients (amplitude along the eigen-PSF), fit 3rd order 
polynomials to the spatial variation of the coefficients, 
and generate an effective PSF model after interpolat- 
ing or stacking the PSF models for the individual ex- 
posures. Twenty principal components or eigen-PSFs 
per exposure can robustly reproduce the observed vari- 
ation of the PSF ellipticity and size (<99% of the vari- 
ance) within each CCD. We measure galaxy semi-major 
and semi-minor axes in the co-add image by fitting el- 
liptical Gaussians convolved with the PSF m odel using 
methodology as in Bernstei n fc Jarvisl ([2002) . Galaxies 
for which this fit converges poorly are not used in the out- 
put catalog of ellipticities. This Stack-Fit measurement 
method has been tested on sim ulations of the Large Syn- 
optic Survey Telescope fLSST; Uee fc Tvsod l20Ll. and 
full details are given in (|Jee et aL^OfHT The difference 
from the procedure described in pee fc Tyson! 1201 fl ) is 
that the images are matched to the DLS: the simulated 
shear maps are convolved with spatially varying PSFs 
similar to what is found in the DLS , down-sampled to 
the DLS pixel scale, and have noise added to match the 
DLS depth. In summary, since the derived shear is in- 
creasingly underestimated with increasing noise, a S/N 
or magnitude-dependent calibration factor is necessary 
to recover the true shear from the measured shear. The 
true shear -y true can be related to the observed shear 7° bs 
using both a multiplicative component, m 7 and an addi- 
tive component C via <-y tmc = TO 7 7° bs +C. Based on the 
above image simulations, the multiplicative calibration 
factor m 7 can be parameterized by: 

m 7 = 6 x 10~ 4 (m fl - 20) 3 26 + 1.036 (11) 

The additive calibration factor C is negligible, as it is 
10% of the statistical error on all radial scales relevant to 
this work. In this analysis, we use source galaxies with 
22<m fi <24.5, corresponding to 1.04<m 7 <1.12. Addi- 
tionally, we require that the semi-minor axis be greater 
than 0.4 pixels to eliminate very small galaxies that 
might suffer from pixellation effects, a e < 0.3, and STA- 
TUS=1, where STATUS is the convergence indicator for 
the minimization ro utine used for fit ting the elliptical 
Gaussians, MPFIT (|Markwardtl I2009D . For our data, 
shape shot noise <7sn ^0.3 (see Equation [7]) . 

Given the position of a foreground lens, the tangen- 
tial and 45°-rotated ellipticity components are calculated 
from the semi-major and semi-minor ellipticity compo- 
nents ei and e2 as follows: 

e t = — ei cos 28 — e2 sin 20 (12) 
e x = ei sin 20 — e2 cos 26 (13) 

where 9 = arctan((y source - yi ens )/ (^source - x\ sns )). 
(xi C ns, yions) and (x sourco , y SOU rcc) are the coordinates of 
the lens and source, respectively. 

The average 45°-rotated ellipticity component is calcu- 
lated as a systematic test as it is expected to be consistent 
with zero at all radial scales. Two additional checks for 
systematics include measuring the lensing signal of stars 
and the lensing signal of source galaxies around random 
locations in the catalog. The results of these systematics 
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checks, all consistent with a null result, are shown along 
with the mean galaxy-galaxy lensing signal for all galaxy 
lenses in the redshift range 0.35 < z < 0.55 and absolute 
magnitude range -22 < Mr < -19 in Figure [7] 

3.2. Photometric Redshifts 

We compute photometric redshifts (photo-z) by fitting 
the four-band photometry to a set of galaxy templates 
to determine both redshift and galaxy type. This is 
done with the publi cly available B ayesian Photometric 
Redshift code (BPZ: lBemtezll2000f ). The code utilizes a 
Bayesian prior for the probability of the redshift given the 
type and magnitude. Instead of using the default BPZ 
magnitude prior which treats galaxies with R<20 iden- 
tically and is based on a small nu mber of galaxies in th e 
Hubble Deep Field North (HDFN: IWilliams et al.lll996h . 
w e fit new param eterized priors following the prescription 
in lBem tcz (2003). For R<21, we calibrate the prior with 
spectroscopically confirmed galaxies in the Smithsonia n 
Hectospec Lensing Survey fSHELS: iGeller et al.ll2005h . 
which is complete to R=20.7. For fainter galaxies, we fit 
to th e VIMOS-VLT Deep Survey (VVDS; lLeFevre et al.l 
2005), which consists of more than 11,000 spectra from 
17.5<iab <24.0. While the possibility remains that there 
is a population of galaxies systematically missing from 
our magnitude complete training sets that could intro- 
duce a bias in the calculated prior, we only broadly fit 
for three galaxy classes (Elliptical, Spiral, and Starburst) 
and use a parameterized N(z), which mitigates small 
missing populations. The resulting prior is qualitatively 
very similar to the HDFN prior but is based on almost 
100 times as many spectroscopic redshifts and has a more 
sophisticated description for R<20 . We empirically ad- 
just the six standa rd CWW+SB (|Coleman et al.l 119801: 
IKinnev et al.l Il996l ) templates using the photometry of 
galaxies with known spectroscopic redshifts (spec-z) in 
SHELS, and we employ the resulting templates to de- 
termine photo-z and K-corrections to z = 0. In Fig- 
ure [T] we show a number density plot of spec-z vs photo- 
z for 10,000 galaxies in the southern field F5, which has 
overlap with the PRIsm MUlti-object Survey (PRIMUS; 
ICoil et al.|[20TTI ). The quoted 100% sampling range is 
R<22.8, the 30% sampling range is 22.8<R<23.3, and 
the redshift precision is a spec - z = 0.005(1 + z trU e-z) due 
to the low spectral resolution of the instrument. The 
root-mean-square cr Zirms of the difference between the ob- 
served redshift and the true redshift of a galaxy at each 
redshift bin is defined as: 



^obs Ztr 



1 



Ztr 



(14) 



After outlier rejection, the DLS photo-z scatter 
o~photo~z = 0.06(1 + z S p ec - z ), and about 4% of the galax- 
ies have photo-z outside of 0.2(1 + z) (catastrophic out- 
liers). Further details can be found in Schmidt & Thor- 
man (2012, in prep.). 

BPZ reports both a redshift probability density func- 
tion (PDF) and a one point summary statistic zb, the 
peak of the PDF. The PDF is the posterior distribution 
given the data, in a Bayesian sense, marginalized over 
template type and apparent magnitude. In this work, 
we use zb when calculating Equation [4] We also use zb 
in order to calculate absolute magnitudes and rest-frame 




0.4 0.6 

PRIMUS spec-z 

Fig. 1. — Density plot of spectroscopic z vs photometric z. The 
photo-z are scattered about the black solid line that represents 
the spec-z=photo-z relation with cr p hoto—z = 0.06(1 + z spec — z ). 
Regions where there is notable bias include galaxies at spec-z~0.1 
that are scattered to larger photo-z and galaxies at 0.3<spec-z<0.5 
that are assigned to photo-z~0.25. This plot shows that we must be 
careful with photo-z<0.3 and take relatively wide vertical photo-z 
slices in our analysis of at least 0.2. 
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Fig. 2. — Photo-z vs absolute R magnitude for redshift and mag- 
nitude selected lenses. The black contours indicate the number 
density, with levels given by logarithmically spaced bins from 1 
to 14,000. The dashed blue lines correspond to the faint apparent 
magnitude cutoff mji max = 24 and the bright apparent magnitude 



cutoff nip , 



18. The black rectangles show the cuts we use in 



Section 1 5 . 3 1 when we make comparisons between redshift shells, la- 
beled zl-3. 



luminosities for the lens galaxies, differentiate the lens 
sample from the source sample, and weight the galaxies 
by the lensing kernel. We use the terms zb and photo-z 
interchangeably for the remainder of the paper. 

Table Q] lists the redshift and magnitude selections ap- 
plied to create lens and source samples. All source galax- 
ies have rxiR < 24.5. This conservative faint cut-off cor- 
responds to the peak of the apparent magnitude num- 
ber counts for the joint catalog containing galaxies with 
both shapes and photo-zs. After further conservative 
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TABLE 1 

Summary of the cuts applied to select the lenses L1-L3, 

Z1-Z3 AND THE CORRESPONDING SOURCES SL AND Szl-Sz3. "-" 
MEANS MINIMUM AND " +" MEANS MAXIMUM, "n": NO THRESHOLD 
IMPOSED. 



cuts in photo-z and magnitude described in Table [T] as 
well as data quality and S/N cuts, the trimmed source 
catalogs used in this analysis contains between 225,580 
and 506,241 galaxies. The corresponding luminosity- 
binned lens sample contains 436,943 galaxies. We plot 
number density contours for photo-z vs absolute R mag- 
nitude in Figure for the redshift range 0.1<z#<0.8 
and 18<niR <24. The three black squares indicate the 
volume-complete lens samples that we later compare in 
Section 15.31 

To check our lens and source selections using zb, we 
show in Figure [3] the arbitrarily normalized summed 
p(z) from BPZ for foreground lenses selected in a fixed 
absolute magnitude range -19<M^<-22 and split into 
redshift shells of 0.15<z B <0.35, 0.35<z B <0.55, and 
0.55<zs<0.75 and for background source samples se- 
lected such that their peak p(z) occurs at more than 0.3 
higher redshift than the peak in the source sample (see 
TableQ]for the specific cuts). The resulting summed p(z) 
show lens distributions that are generally well-separated 
from the source distribution. We consider the effects of 
photo-z inaccuracy including the residual systematic er- 
rors due to the ~4% outliers on the lensing analysis in 
Section O 

3.3. Star/Galaxy Separation 

For star/galaxy separation, we apply two cuts. At 
bright magnitudes, stars are easily distinguishable from 
galaxies in a size vs magnitude diagram. We use SExtrac- 
tor FLUX_RADIUS as our measure of size, and cut out 
objects with FLUX_RADIUS<2.4 pixels with niR <19. 
The pixel size is 0.257 arcseconds. At fainter magni- 
tudes, the star and galaxy loci overlap on the size vs 
magnitude diagram, so we employ a cut which involves 
all three second central moments of the PSF and the 
object. For each object, we determine the total % 2 for 
the hypothesis that it is a point source, using the three 
second moments of the object and of the PSF at that lo- 
cation. We then cut on Xtot > 5 to eliminate most faint 
stars. This procedure was shown to result in ~ 0.4% stel- 
lar contamination in spectroscopic targets for m# < 22.5 
([Dawson et al.ll2~011[ L At fainter magnitudes, the num- 
ber density of galaxies overwhelms that of stars while the 
efficiency of the \ 2 cu t should not diminish (although 
a larger fraction of galaxie s might be cut). U sing stel- 
lar population models from Robin et al. (2003) and DLS 
coordinates to estimate star counts and typical galaxy 
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Fig. 3. — The summed normalized likelihood function p(z) from 
BPZ for redshift and magnitude selected sets of lenses and sources. 
The cuts applied to create the lens and source samples are de- 
scribed in Table [T] The lens and source samples are generally well- 
separated. However, there are non-negligible overlaps between the 
z2 lenses and sources and between the zl lenses and z2 lenses. We 
address the effect of these overlaps on our analysis using mock cat- 
alogs, and angular cross-correlation s wi thin the observational data, 
and discuss our findings in Section |4. II 

counts from the literatureQ we estimate the stellar con- 
tamination to be < 0.2% for mji ~ 24. The shape cuts 
described in Section 13.11 further reduce this stellar con- 
tamination. 

4. PHOTOMETRIC REDSHIFT ERROR ANALYSIS 

Photo-z errors can affect the amplitude of the galaxy- 
mass correlation in several ways, which we can separate 
into those related to errors in lens galaxy photo-z (zl) 
and those related to errors in source galaxy photo-z (zs)- 
First consider errors in lens photo-z. Error in zr, will 
smear the signal in projected comoving radius due to 
error in the conversion from an angular distance to a co- 
moving distance. Especially for sources close to the lens 
redshift there will also be an error in AE due to an er- 
ror in E cr j t (see Equation [7]) . In the case that zl > z$, 
the signal will be diluted by "source" galaxies which are 
in reality foreground objects. There is an additional is- 
sue that the errors in zl propaga te to errors in absolut e 
magnitude Mr. As found by e.g. iHoekstra et a l. (2005), 
the net effect (due to the shape of the galaxy number 
counts) is that intrinsically fainter galaxies are more of- 
ten upscattered into brighter luminosity bins, causing a 
dilution in the observed signal. Now consider errors in 
source photo-z. For a source galaxy that is behind a lens 
galaxy an error on a source redshift will create a bias 
in AE due to an error on E cr i t . An error in a source 
redshift such that zs < zl will create a signal dilution 
because this "source" galaxy is now a foreground object. 
The net correction is a mix of all these effects, and the 
mix changes with galaxy a nd lens samples. For fu rther 
discussion s ee Figure 8 of | Leauthaud et al.l (|2010l ) and 
Section 7 of lNakajima et all (|2012t l. 

We undertake a galaxy-galaxy lensing Monte Carlo 
simulation, creating mock lens and source catalogs rep- 
resenting DLS-sized fields of 4 deg 2 to understand 

2 See |http:// star- www, dur .ac.uk/ ~im /pub ht ml /counts /counts . html 
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how photo-z scatter, photo-z catastrophic outliers, and 
lens/source sample selection affect the accuracy of the 
projected galaxy-mass correlation measurement. We also 
compute angular cross-correlation between the lens and 
source samples to investigate the purity of these samples. 

4.1. Signal Recovery Simulations 

We calculate the ratio of observed (with scat- 
tered redshifts) lensing signal to input lensing signal 
A£(i?) scat t/AE(i?)j n p Ut for each of the three lens sam- 
ple redshift cases. For the mock source galaxies we draw 
random x and y positions for simulated source galax- 
ies, and we draw redshifts and intrinsic ellipticities in a 
Monte Carlo fashion from distributions consistent with 
the data. We draw the intrinsic ellipticity components 
along the major and minor axes from Gaussians with 
width 0.3 and centered at 0. The summed observed p(z) 
for sources (as shown in Figure [3]) is a broadened ap- 
proximation for the true dN / dz of the sources since this 
distribution already includes observational effects such 
as photo-z errors and is broader than the underlying dis- 
tribution. To simulate the native distribution of galaxies 
in redshift for each galaxy population, we find that dis- 
tribution p0(z) which when convolved with the known 
o~ z photo-z errors gives the observed summed p(z) for 
that population. We draw redshifts from the summed 
intrinsic p0(z) distribution of the sources. 

For the mock lens galaxies, we fix the x, y, and luminos- 
ity of the simulated lens galaxies to the values measured 
from the observed catalogs. That is, each simulated lens 
galaxy has a one-to-one counterpart in the real catalog 
(e.g. that of F5) from which it inherits x and y coordi- 
nates in pixels and luminosity. These mock lens galaxies 
thus inherit the angular clustering that is present in the 
observed sample of lens galaxies. As with the sources, 
we draw redshifts for the lenses from an estimate of the 
intrinsic summed p0(z) of the lenses which when con- 
volved with the photo-z error yields the observed p(z) as 
shown in Figure [3l We also take into account any red- 
shift overlap between the lenses and sources and among 
the zl, z2, and z3 lens bins. To convert the lens lumi- 
nosities to halo masses consistent with the data, we fit 
NFW profiles to the inner 300 kpc of the lensing signal 
in each luminosity bin and fit a power law to the mean 
luminosity vs halo mass as shown below in Figure [5j The 
best fit mass-luminosity relation is: 

Af 2 oo = 1-46 x 10 3 Af Q (L/L G )°- S9 (15) 

To mitigate edge effects, we pad the boundaries of each 
simulated mass field with an additional degree on each 
side. That is, the lens galaxies are distributed over a 
total area of 16 deg 2 for each of the five fields. 

Using NFW mass profiles for the lens halos, we cal- 
culate the expected shear for each lens-source pair given 
the physical separation and lens Af2oo- For each source 
galaxy i, we use the weak lens approximation in calcu- 
lating the total applied shear due to all lens galaxies 
7T,i — X^ lona 7j • The observed ellipticity is then given 
as in Equation [3] The resulting mock lens and source 
catalogs have the selection effects of our observations, 
and is a starting point from which to test how photo- 
z errors such as scatter, photo-z bias, and catastrophic 
outliers affect the galaxy-mass correlations reconstructed 



from the lensing analysis via the corresponding errors 
in the distance ratios. To examine this, we assign "ob- 
served redshifts" by applying photo-z errors similar to 
those measured from the observational data. We obtain 
the observed redshifts by drawing from a Gaussian cen- 
tered at the assigned true redshift with width given by 
a value describing the photo-z scatter (see Section l3~2"j) . 
We use a photo-z scatter of cr z , rms =0.06(l+z). We also 
introduce catastrophic outliers by switching the redshifts 
of a random 4% of the galaxies with random redshifts. 

After determining the absolute magnitudes corre- 
sponding to the scattered redshifts, we re-select lenses 
and sources using these new "observed" redshifts and ab- 
solute magnitudes and measure A£(i?) sca tt C r as in Equa- 
tion!?] At very large angular separation there is no mech- 
anism that can cause a spatial dependence of the ratio 
A£ (-R)scatt/ AS (R) input) and we fit a constant over 1-10 
Mpc radius. We repeat this entire simulation for various 
levels of photo-z error and find that the recovery ratio 
for each lens redshift sample decreases with photo-z er- 
ror. The best fits for this recovery ratio for 0.06(l+z) 
scatter are 0.81, 0.85, and 0.90 for LI, L2, and L3, re- 
spectively and 0.94, 0.83, and 0.74 for zl, z2, and z3, 
respectively. In the analysis in Section 7.3 we correct for 
this systematic error. 

4.2. Angular Cross-Correlations 

Do the summed ~p(z) distributions for our lens and 
source samples represent the true distributions? While 
the distributions of PRIMUS spec-z, DLS photo-z, and 
p(z) are consistent out to the limit of the spectroscopy at 
z ~ 1, it is important to investigate consistency for our 
entire source and lens sample galaxies in other ways. As 
a test of source and lens sample purity we calculate angu- 
lar cross-correlations wi 2 (6 l ) between these samples. Wc 
find that the ratios between the auto and cross correla- 
tions of these samples is less than 0.1, consistent with the 
overlap of the p^i^) tails and the known photo-z outlier 
rates. The amplitude of wi2(#) never rises above 0.01. 
The residual sample impurity results in a slight decrease 
in signal-to-noise ratio, but is small enough that we do 
not make any additional corrections in this analysis. 

5. RESULTS 

We measure the cross-correlation between the lens 
galaxies and the associated total projected mass distri- 
bution which is in the foreground of the source galaxies. 
The DM halo in which the galaxy is resident, as well 
as all other projected mass correlated with the galaxy, 
contribute to this weak lens signal. We bin the shear sig- 
nal in logarithmic bins by projected comoving radii and 
convert to a projected differential surface mass density 
following Equation [71 Errors are calculated via jackknifc 
resampling where the jackknife components are one ninth 
of a field size. Over the 5 fields, this yields a total of 
45 jackknife subsamples. These data permit an analysis 
of the distribution of mass associated with foreground 
galaxies over a wide range of projected physical scales, 
and its correlation with galaxy luminosity and redshift. 

5.1. Galaxy-Mass Correlation vs Luminosity 

We divide the lenses into three ranges of absolute mag- 
nitude for a fixed range of 0.35<zs<0.75 and measure 
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Fig. 4. — Galaxy-mass correlation (expressed as a mass over- 
density around lens galaxies) vs radius and luminosity. The sur- 
face mass overdensity for three bins of lens absolute luminosity 
as a function of radial separation is shown. All galaxies have 
0.35<zg<0.75. Red squares show Ll, green triangles show L2, 
blue circles show L3, and the NFW halo model fits to the inner 
300 kpc are overplotted as dashed lin es. C orrection factors from 
photo-z errors, as discussed in Section 14.11 have been applied. L2 
and L3 data points are slightly horizontally offset for clarity. On 
small radial scales, the profiles scale with luminosity such that the 
brighter lensing galaxies have higher galaxy-mass correlation. On 
large radial scales, the profiles flatten, indicating additional mass 
contributions from neighboring halos. The shape of these radial 
profiles and trends with lumin osity are similar to the t heoretical 
prediction shown in Figure 8 of Hayashi & White (2008). 



the galaxy-galaxy lensing signal for each lens galaxy lu- 
minosity sample. The properties of these samples are 
summarized in Table [TJ After dividing the signal by the 
correction factors due to photo-z errors discussed in Sec- 
tionO(0.81, 0.85, and 0.90 for Ll, L2, and L3, respec- 
tively), the results are shown in Figure 2] along with NFW 
profiles fit to the inner 300 kpc of each A£(i2). These 
fits take into account the full covariance matrices esti- 
mated from the data (see the Appendix for details about 
how the covariance matrix is calculated), although this 
inclusion does not have a significant effect due to our 
focus on small radial scales. The signals scale with av- 
erage luminosity with higher surface mass overdensities 
corresponding to brighter average luminosities. As ex- 
pected, intrinsically more luminous galaxies have higher 
mass profiles inside their virial radii than low luminosity 
galaxies. 

As shown in Figure 0] on scales larger than ^0.3 Mpc 
all three signals flatten, and converge on 10 Mpc scales 
to a narrow range of mass overdensity. This is because at 
large radii the galaxy-mass correlation reflects the mass 
auto-correlation modulo t he galaxy bias, as d i scuss ed in 
IHavashi fc Whitel (l200l) . IHavashi fc Whltel (120081) use 
the Millenium simulation to examine the cross correla- 
tion between halo centers and mass, showing a depen- 
dence on galaxy mass (or luminosity) similar in profile 
shape to our results (their Figure 8). The strong lumi- 
nosity dependence of AS at R<1 Mpc and weak lumi- 
nosity dependence at R>lMpc have also been found in 
the simulati ons of Tasitsio mi et al.l ()2004l ) and the ob- 
servations of lSheldon et al.l ( 20041 ). The behavior of the 
galaxy-mass correla tion at large r adii h as also been re- 
cently examined bv lMasaki et al.l (|2012l ) in a 1024 3 par- 



P(m ha i \l) oc exp 



where I 



ticle N-body CDM simulation. 

5.2. Mass vs Luminosity 

From the best fit NFW profiles, we obtain values for 
M200 for each luminosity bin. It is worth noting that the 
NFW is not a good fit to the AE profiles beyond a few 
hundred kpc. In particular, the L3 sample likely con- 
tains a non-negligible fraction of satellite galaxies, and 
preliminary work fitting a full halo model indicates M200 
for this luminosity bin is overestimated by ~20%. 

As pointed out in iTasitsiomi et ail (|2004D . the best fit 
NFW mass generally falls between the mean and median 
mass for a broad halo mass distribution. The conversion 
from the best fit NFW mass to the mean mass is an up- 
wards correction that depends on the scatter in the mass- 
luminosity relation, which is higher at the bright end 
where the halo mass distribution is broader. To deter- 
m ine this correction, we follow the procedure described 
in Ivan Uitert et ail ()201lD . We adopt a conditional prob- 
ability function for the halo mass given a luminosity of 
the form: 

/ (™halo — ™halo,cent) 2 \ 
I J (16) 

log(L), TOhaio = log(M ha i ), and cr^ ha]o is the 
scatter in mhaio- For a given best fit NFW mass and a 
(Tjjjj, from work on satellite kinematics in SDSS bv lMorel 
(2011) (their Fi gure 4 produces values from 0.35 to 0.45 
for our mean luminosities), we then convolve the condi- 
tional probability f unction with the mass function from 
iTinker et al] pOOl and draw 1000 masses from the re- 
sulting distribution. We calculate and average the NFW 
A£(i?) profiles for this ensemble of masses and compare 
the best fit NFW M 2 oo to (M 20 o}. We obtain conversion 
factors of 3.8 (Ll), 3.1 (L2), and 2.5 (L3) and multi- 
ply our best fit NFW M200 by these numbers to obtain 
(M 20 o)- 

We plot (A/200) vs the mean luminosity of each bin in 
Figure [5] along with a simple power law fit (Equation ll5[) . 
In Figure El our mean luminosities are converted to the 
AB system. Since the mean redshifts of the three lumi- 
nosity bins decrease slightly with decreasing luminosity, 
we also scale the DLS masses to the mean redshift of 
Ll (z=0.59) for the sake of consistency. The mean red- 
shifts of L2 and L3 are 0.58 and 0.51, which corresponds 
to multiplication factors of 1.01 and 1.10 based on our 
M200 definition which depends on redshift through p cr it- 
We perform the power law fit primarily for use in the 
Monte Carlo simulations described in Section 14.11 We 
also show results for early-ty pe (ET) and late-type (LT) 
gala xies at low redsh i ft fro m IMandelbaum et al.l (2006) 
and Ivan Uitert et al.l (|2011l ). However, there are several 
caveats to consider befor e making a comparison. First, 
IMandelbaum et all (|2006l ) use a different mass definition 
(180/5 instead of our 200/9 c ), so our masses should be 
adjusted upwards by ^30% f or a more direct compari - 
son. In Figure we plot the IMandelbaum et all (|2006f ) 
masses after dividing by 1.3 to make the qualitative com- 
parison. Second, the DLS points are at a redshift of 0.59, 
and we have not yet included corrections for passive evo- 
lution. If we adopt a passive evolution correction from 
iBlanton et al.l (|2003l ) of 1.6(2; — 0.1) for galaxies that are 
best fit with elliptical templates (e.g. BPZ T_B<1.5), 
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Fig. 5. — The mean luminosity vs mean halo mass for the three 
luminosity bins described in Table [T] Halo masses calculated by 
fitting the NFW profile to the inner 300 kpc are shown as bl ack 
stars and then using the conversion factors described in Section l5,2l 
to obtain mean masses. The dashed red line shows a pow er law fit 
to the mean masses. The power law given by Eq uatio n lf 5l is used in 
the Monte Carlo simulations discussed in Section [4.1l to convert lens 
luminosities to halo masses. For comparison purposes, previous 
results ar e also plotted, although no te that different methodologies 
are used. Mandelbaum ct al. (2006) results at (z) ~0.f are shown 
as blue triangles for late-types (IT) and thin magenta diamonds 
for e arly-types (ET). Note tha t we have divided the data points 
from[Mandclbaum ct al. (2006) by 1.3 to correct for different mass 
definitions. Ivan Uitert et al.l 1120111 ) results at (z) ~0.1 are shown 
as cyan filled circles for late-types and yellow squares for early- 
types. The errors on the SDSS results are 2-cr. The DLS results 
agree well with those in the literature after making adjustments 
for mass definitions and passive luminosity evolution as discussed 
in the text. 

K-correct to z = 0.1, a nd calculate the dis t ance m odu- 
lus with h = 1.0 as in Ma ndelbaum et all (|2005aD . our 
mean luminosities decrease by ^30%. Third, we have 
not split our lens samples by type, so the DLS masses 
include both ET and LT galaxies. Finally, the conver- 
sion factors we applied to obtain the conversions from 
best fit NF W Mo p n to ( Af 2 oo) are based on our adopted 
cr mhalo from lMorel (|2011l) whose results are limited to z~0 
central galaxies. The true c mhalo for our higher redshift 
central and satellite galaxies might be somewhat differ- 
ent, ultimately changing our interpretation of (M2oo)- 

5.3. Galaxy-Mass Correlation vs Redshift 

We next divide the lenses into three ranges of redshift 
for a fixed range of absolute magnitude -22<Mr <-19 
and measure the galaxy-galaxy lensing signal for each 
redshift sample. This range of absolute magnitude was 
chosen such that the samples could be as volume com- 
plete as possible (see Figure [2J while still allowing a rea- 
sonable S/N. The properties of these samples are sum- 
marized in Table [TJ and the results are shown in Figure H] 
after applying the corrections for photo-z errors found in 
Section|LTJ 0.94 (zl), 0.83 (z2), and 0.74 (z3). The bins 
at large separations are correlated, and we discuss how 
we calculate the covariance matrices in the Appendix and 
show the normalized covariance matrices in Figure [9j At 
radii less than a few hundred kpc, all three signals have 
similar amplitudes, indicating that the mass range corre- 
sponding to the luminosity selection is consistent across 
redshift shells. However, there is a noticeable trend at 




R (Mpc) 

Fig. 6. — Redshift dependence on all scales for a single wide 
range of luminosities -22<Mj; <-19. We compare the galaxy-mass 
correlation vs radius for three lens redshift shells, which have prop- 
erties described in Table [l] These signals include corrections for 
photo-z e rrors that were derived from mock catalogs as discussed 
in Section l4.ll zl is shown in black, z2 in magenta, and z3 in cyan. 
The z2 and z3 points are slightly horizontally offset for clarity. On 
small radial scales, the three samples reflect signals consistent with 
a similar mass. However, on larger scales, there is a trend where 
higher AE corresponds to lower mean redshift. 

larger radii where the amplitude of AE increases with 
decreasing redshift. Both before and after applying the 
recovery corrections from Section 14.11 there is a trend. 

6. DISCUSSION 

We have presented projected galaxy-mass correlations 
measured using galaxy-galaxy lensing over a wide base- 
line in radius, luminosity, and redshift. On small radial 
scales, these signals probe the individual DM halos in 
which galaxies reside. On larger radial scales, these sig- 
nals are sensitive to correlated large-scale structure and 
neighboring DM halos. We have focused on two lines of 
investigation: fixing the redshift range and varying the 
mean rest-frame luminosity and conversely, fixing the lu- 
minosity range to a volume-complete sample and varying 
the mean redshift. 

For fixed redshift and varied luminosity, we find the 
well-established scaling of the halo mass with luminos- 
ity: intrinsically brighter galaxies are also more mas- 
sive. In Figure [SJ we show a comparison of our re- 
sults for the M^ a ; - L rel ation at z~0.6 with t h ose o f 
iMandelbaum et al.l (|2006D and Ivan Uitert et al.l (|201lD 
at z~0.1. After corrections for differing mass defini- 
tions and passive luminosity evolution, we find reason- 
able agreement with the literature results. 

On larger radial scales, the shape of the galaxy- 
mass correlations and their dependence on galaxy lu- 
minosity look qui t e sim ilar to the simulation results of 
lHavashi fc Wh ite (2008). There are currently few lens- 
ing observations of galaxies over wide ranges of luminos- 
ity and redshift at radii greater than a few Mpc. We can 
qualitatively compare our results shown in Figure U with 
four observations at lower redshift. Usi ng shear mea- 
surem ents around z~0.1 SDSS galaxies, iSheldon et al.l 
(2004) found trends with luminosity that are similar to 
the higher redshift findings presented here; however, they 
fit power laws to their shear and galaxy-mass correlations 
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and did not see significant deviation from power laws 
at large sca l es. W e can also compare qualitatively to 
iReves et al.l (|2010h . who measure shear around luminous 
red galaxies (LRG) in SDSS out to large radii and model 
it with a halo model. Their LRG mass overdensity at 
z=0.3 is consist ent with our low redshi ft zl and z2 sam- 
ples at 10 Mpc. Ivan Uitert et al.l (|201lD study galaxies in 
RCS2 with spectroscopic redshifts from SDSS, limiting 
the lens redshifts to z^O.l. Their Figure 8 shows halo 
model fits to the lensing signal out to 10 Mpc, which 
is consistent with the signal we measure. In these pre- 
vious studies, the lens galaxies are closest in mass to 
our LI sample, but are ge nerally more massive. Finally, 
iMandelbaum et ail (|2012f ) show results for AS for three 
redshift bins out to z~0.5. However, a direct comparison 
is difficult because they probe much more massive and 
more highly biased galaxies (i.e. LRGs). 

Most interestingly, our galaxy-galaxy lensing data al- 
lows an investigation of the redshift evolution of the pro- 
jected galaxy-mass correlation on 10 Mpc scales at fixed 
lens galaxy luminosity. We find evidence for growth over 
time of the galaxy-mass correlations on large scales from 
z=0.65 to z=0.2 as shown in Figure [6] Interpretation of 
this result in terms of evolution of bias and growth of 
LSS mass structure with cosmic time requires additional 
information. The overall effect as a function of lens sam- 
ple redshift depends on both LSS growth and galaxy bias 
as a function of mass and type, and an additional probe 



such as the lens galaxy clustering signal is necessary to 
disentangle the growth and bias. We will combine this 
two-point auto and cross-correlation information for each 
of the three lens samples with the galaxy-mass correla- 
tions reported here in a joint analysis in future work. 
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Fig. 7. — The mean tangential shear for the z2 lenses and three systematics checks as a function of angular scale. Note that the y-axis is 
log-scale down to 0.005 and linear-scale below in order to magnify the small systematics test results. The black circles represent the signal 
from 45°-rotated source shears, the red circles show the signal from one million random positions cross-correlated with star shears, and 
the green circles show the signal from one million random positions cross-correlated with source galaxy shears. The data points are offset 
horizontally for clarity. Errors are estimated from jackknife resampling, and the errors on the random-star cross-correlations are too small 
to be visible on the plot. The systematics tests are all consistent with zero, and generally at least an order of magnitude below the mean 
tangential shear of the real foreground galaxy positions cross-correlated with background galaxy shears (magenta circles). 



APPENDIX 
SYSTEMATICS TESTS 

We undertake several tests for shear and photo-z systematics. Figure [7] shows the results of three checks for shear 
systematics as a function of angular scale: foreground lens positions cross-correlated with tangential shear from 
45-degree rotated source galaxies, one million random foreground positions (200k per field) cross-correlated with the 
tangential shear from source galaxies, and one million random foreground positions cross-correlated with the tangential 
shear from stars. In order to plot the positive and large shear signal on the same plot with the relatively small residual 
bipolar systematics, the y-axis is log above 0.005 shear and linear below. These tests are consistent with zero residual 
shear systematics and are typically less than the true signal by at least an order of magnitude. The error bars are 
given by jackknifmg 9 subsamples for each of the 5 fields. The two tests using actual source galaxy shears have the 
largest errors since there is a much larger variation in the galaxy ellipticity distributions than for the stars, which 
should have minimal variation by construction. Note that on large scales, the errors become correlated. 

If cither the photo-z of the lens or source galaxy samples are systematically in error, then the error in the distance 
ratio propagates to an error in the galaxy-mass correlation. Such an error would be revealed by differences in amplitude 
in the galaxy-mass correlations computed using different source samples for the same lens sample. Figure [5] shows the 
cross-correlation of the zl lens sample with four source samples with differing mean redshifts. We find no evidence of 
a distance ratio scaling inconsistency at the level of the noise. 

COVARIANCE MATRICES 

While the systematics in the data have been addressed and reduced to below the noise level, fits of models to the 
data must take account of covariances: the errors can be correlated. The errors in the data shown in Figures |4] 
and [6] are correlated (particularly at large scales) due to common sources used in different lens-source pairs that are 
stacked together. Thus the 13 radial bins are not completely independent in the sense that their errors are correlated. 
It is informative to calculate correlation matrices by re-sampling galaxy-galaxy lensing data for each lens redshift 
population. The correlation matrix is the normalized covariance matrix, with matrix elements given by Corr^ = 
Cij /(Ci^Cj.j) 1 ^ 2 . The covariance matrix C for N jackknife sub-samples is estimated using: 

Y , N 

C (AS (iJi), AS(i£,-)) = -^£[AS fe 0RO - AS^)] [AS* (i?,-) - AS(i?,)] (Bl) 

k=l 

where AS(i£,) = £fe=i AT, k (Ri)/N is the mean value over N sub-samples. For each of the three lens redshift samples, 
we take N = 9 random re-samples of all the galaxies in each field. The correlation matrices corresponding to the 
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Fig. 8. — Lensing of zl by source samples with different mean redshift. We show a test of photo-z and lensing consistency by measuring the 
galaxy-mass correlation in one lens galaxy sample by using all the possible source samples. AS(R) for the zl lens sample is cross-correlated 
against four source shells with very different mean redshifts, described in the legend. The black circles show the cross-correlation with the 
zl background sources as shown in Figure 7, the red circles show the cross-correlation with z3 lenses (used in this test as a background 
sample), the cyan circles show cross-correlation with the z2 background sources, and the yellow circles show the cross-correlation with the 
z3 background sources. Table 1 describes the cuts for each of these samples. The data points are calculated for the same radial bins but are 
offset for clarity. While the four signals are correlated with each other, the mean redshift of the background sample varies. The distance 
ratio scaling appears to be consistent. 



lensing signals for the three redshift bins in Figure [6] are shown in Figure [9l Each cell represents the level of correlation 
between the given pair of radial bins with darker shades corresponding to higher correlation. There is a higher level 
of correlation at large radial scales at low lens redshift. This may be caused by the correlated noise in the dithered 
40' subfield observing. At higher redshift 10 Mpc is completely within the subfield angular scale. Scales of 1 Mpc 
correspond to angular scales of 4.3' for lens sample zl, 2.9' for z2, and 2.4' for z3. At the higher redshifts of z2 and z3, 
there is an increasing number of lens-source pairs within a given 40' subfield due to the conversion between angular 
and physical scales. However, note that due to our cuts the overall number density of source galaxies decreases going 
to higher lens redshifts (refer to Table[T]for exact numbers). Thus, the stronger correlations for zl are also likely linked 
to the fact that zl has the largest number density of source galaxies. While the correlation between bins does not 
affect the actual values of the galaxy-mass correlation data points, the correlated errors need to be taken into account 
in fits. We consider the full correlation matrices in the NFW fits presented in this work. 
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Fig. 9. — Correlation matrices for the lensing signals from the three redshift shells shown in Figure[6] The leftmost panel for zl reflects 
higher levels of covariance going to large radial scales, zl also has the largest number density of source galaxies, and the overall source 
number density decreases for z2 and z3. Going to higher redshifts, a given angular scale corresponds to a larger physical scale and the 
covariance at large radius moves to the upper right. At higher redshifts, there are more lens-source pairs within a given subfield due to the 
angular to physical scale conversion, and the subfield angular scale moves to projected scales much larger than 10 Mpc. The off-diagonal 
cells in the panel corresponding to z3 are mainly noise, since most of the pairs come from within a subfield. 



